Condensation and vortex formation in Bose-gas upon cooling 
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The mechanism for the transition of a Bose gas to the superfluid state via thermal fluctuations is 
considered. It is shown that in the process of external cooling some critical fluctuations (instantons) 
are formed above the critical temperature. The probability of the instanton formation is calculated 
in the three and two-dimensional cases. It is found that this probability increases as the system 
approaches the transition temperature. It is shown that the evolution of an individual instanton is 
impossible without the formation of vortices in its superfluid part. 

PACS numbers: 



I. INTRODUCTION 

The ideas of the kinetics of phase transitions have 
been thoroughly developed for first-order phase transi- 
tions and envisage the existence of the metastable phase 
itself and an equilibrium critical nucleus. The corre- 
sponding theory was worked out in , 0] and described 
in detail in However, theoretical concepts concern- 
ing the kinetics of second-order phase transitions, where 
these two facts do not exist, have been developed insuf- 
ficiently. In the was proposed a certain special model 
for the formation of an ordered phase after the fast phase- 
transition stage in the " short-range" order in the presence 
of only two types of ordering. 

The interest in the problem of a phase transition upon 
a fast change in external parameters (e.g., temperature) 
has been aroused in connection with the cosmological 
ideas of the Big Bang, where the rapidly expanding Uni- 
verse must be cooled and pass through a series of phase 
transformations accompanied by a change in the symme- 
try of physical fields 5] . It was proposed that the kinet- 
ics of these transformations can be modeled in condensed 
matter |6j. 

In the '?] was proposed a theory of the second-order 
phase transition upon a rapid change in temperature in 
liquid 4 He. The main assumption in the proposed mech- 
anism is about the "critical retardation" of all processes 
in the vicinity of the transition temperature and " fast" 
formation of the nuclei of a new phase upon the subse- 
quent cooling. This gives rise to a large number of defects 
on the order of the number of fluctuations far above the 
transition point. 

However, no retardation in the formation of a new 
phase has been detected experimentally; the critical re- 
tardation is associated with the duration of the equilibra- 
tion process at macroscopic distances, which is insignif- 
icant for the nonuniform process of formation of a new 
phase. 

In this work, we consider the transition to a new phase 
via the evolution of fluctuations on scales much smaller 
than the correlation length, which can occur quite rapidly 
even in the vicinity of the critical temperature. The tran- 
sition kinetics in this case is found to be directly related 



to the cooling process itself. In our preliminary paper [8j 
this approach was proposed for the specific problem of 
Bose condensation of a weakly interacting Bose gas. The 
appropriate set of equations governing critical fluctua- 
tions was derived. In the present work a detailed analysis 
of the main equations is performed including numerical 
calculations and a generalization to the case of the two- 
dimensional exciton gas. We will consider the formation 
of a condensate in the model of a weakly nonideal Bose 
gas with external cooling and demonstrate an analogy 
with first-order phase transitions. 

A similar approach to the problem of the wave nucle- 
ation rate due to thermal noise was considered in [jj. 
But the problem of the present work requires the noise 
to be connected to the random thermal fluxes. In this 
case it is appropriate to use a more general approach for 
the instanton formation namely local Hamilton equations 
rather than the Lagrangian equations used in (see also 
[13). In [HI the problem of large negative gradients in 
Burgers turbulence was considered, which has some re- 
semblance to the differential equations discussed in our 
work, but the boundary problem is quite different. 



II. DILUTE BOSE-GAS UPON COOLING 

The standard theory of a weakly nonideal Bose gas 
involves a Hamiltonian of the form : 



a P3 a P2 a Pl , (1) 



where do - is the scattering amplitude having the atomic 
scale and m - is the atomic mass. The properties of 
such gas for a small density n (determined by the gas 
parameter r\ — naa 3 <C 1 ) are close to the properties of 
an ideal Bose gas with the transition temperature [T^ | : 



T c = 



3.31 h 2 n 2 ^ 3 

y/2 m 



(2) 



At temperatures below the transition point, the ideal 
Bose gas has a pressure depending only on the tempera- 
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ture: 



P ld = 0.0851- 



l 3/2 T 5/2 



(3) 



which corresponds to zero isothermal sound velocity. 

Considering the finite scattering amplitude we can 
write the qualitative equation of state below the tran- 
sition point as 



P = Pd(T) 



h 2 aon 2 



(4) 



We omitted the insignificant constant factor in the second 
term. 

The entire kinetics is essentially determined by the 
Bose-gas cooling mechanism. We will consider a sim- 
ple model where the Bose gas is in a certain solid matrix 
with which it only slightly interacts. Such a situation 
may take place, for example, for the exciton gas in a 
crystal. The crystal can be rapidly cooled to a low tem- 
perature; in this case, the Bose-gas cooling proceeds via 
phonon emission. Assuming that the heat capacity of 
the crystal is large compared to the Bose gas, we can 
disregard the presence of thermal phonons in the crystal 
and their effect on the Bose gas. As a result, we obtain 
a uniform energy-loss mechanism, which is described by 
a phenomenological quantity T/r p h- The other models 
of cooling necessitate the analysis of heat transfer at the 
sample boundaries, which is a much more complicated 
problem. The loss rate l/r p h is determined by the colli- 
sions of particles with each other and by the interaction 
with phonons, which will be regarded as weak: 

l/r p h < l/r tr . 

Since l/r p h ~ nvT<J p n > l/it r ~ nvrao 2 { v t s the thermal 
velocity), this means that 

2 

<y p h *C a , 

where a p h is the cross section for scattering with phonon 
emission, which corresponds to the weak interaction of 
Bose gas with the crystal. 

In view of the smallness of quantity 1 /r p h the evolution 
of the Bose system is slow; in particular, we assume that 
the acoustic wavelength cr p h ~ vxT p h is large compared 
to the characteristic length ~ JX T ph , where x - is the 
thermal diffusivity: 



(5) 




VTTph 



(I stands for the mean free path). This makes it possi- 
ble to assume that the fluctuation evolution occurs at a 
constant pressure that coincides with the thermodynam- 
ically equilibrium pressure. 

It follows from Eq.Q that the density variation Sn in 
the fluctuation region is related to a change in tempera- 
ture by 



Sn 
n 



ST 1 

T fjs 



(6) 



The relative density fluctuation is large compared to 
the relative temperature fluctuation in the temperature 
range T < T c . This leads to a rapid increase in the re- 
ciprocal phonon time 



. 1 ST I 

o — ~ — - 

Tph 



1 



T 



r]3 7p/i 



(7) 



(where l/r p ^° = nvxc p h) with decreasing temperature 
and enhancement of cooling in the fluctuation region. For 
this reason, we will disregard the phonon emission in the 
region far from the developed fluctuations, assuming that 



ST 1 1 



T p h 



T, 



c r]3 T p h 



(8) 



where U(T C -T) = 1 for ST = T-T c < and U(T C -T) = 
for T - T c > 0. 

This allows us to consider the problem of fluctuation 
kinetics within the framework of the theory of hydrody- 
namic fluctuations by supplementing the hydrodynamic 
equations with the energy flux carried away as a result 
of phonon emission: 



T Sn 
r p h° n 



T c -T 

Tph 



U(T C ~T), 



(9) 



where l/T„h = (l/r p ft )(l/?7 1 / 3 ). In view of the con- 
stancy of pressure, we can describe the evolution of tem- 
perature fluctuations by the heat conduction equation 



dT dT 



T-T c 

Tph 



nc p U(T-T), (10) 



where the energy flux carried away by phonons is added, 
x — is the heat conductivity, c p - is the specific heat 
per particle under a constant pressure. This equation 
contains the drift term with mass velocity v(r), which 
appears due to the high density in the fluctuation core. 
In the following analysis, this term will be omitted as a 
higher-order term in fluctuation. We are interested in 
the temperature-field fluctuations and their time evolu- 
tion. To analyze these fluctuations, we must introduce 
random heat fluxes q 0, 01, i.e., the Langevin term 
Vq. These fluxes are delta-correlated (i.e., correlated at 
distances and time intervals smaller than the hydrody- 
namic scales). In the case considered, this is ensured by 
the fact that the time T p h and distance ^fxjph (x ~ is 
thermal diffusivity) are larger than the microscopic char- 
acteristics. 

The probability Wt (T(r)) of realizing the given config- 
uration T{r) fluctuation field at time t obeys the Fokker 
Planck equation in variational derivatives 



ys tr 
.14] 



dt 



ST(r) 



oo y2 



XV 2 T + U(T C - T) 



8T{r) 
T-T,~ 



Tph 



w 



w 



d 3 r. (11) 
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In the absence of the interaction with phonons, the sta- 
tionary solution to this equation coincides with the result 
obtained in the thermodynamic theory of fluctuations. 
The quantity 



X V 2 T + U(T C - T) 



T-T c 

Tph 



dT 

tt- ( 12 ) 
dt y ' 



is the temperature- variation rate upon the deviation from 
the mean value T = . 

We assume that fluctuations occur at a fixed temper- 
ature > T c . Fluctuations with AT = T - T m < 
occur quite frequently and are characterized by a certain 
(in fact, stationary) spatial distribution that determines 
the value of W t (T). In view of the normalization, the 
latter quantity gives the number of small fluctuations in 
a unit volume. However, rare large-amplitude fluctua- 
tions with T^T C — Too, T < T c , also sometimes occur, 
initiating the effective cooling by phonons, so that the 
fluctuation becomes irreversible and the nucleus of a new 
phase appears. Our goal is to calculate the probability of 
such fluctuations in a unit volume per unit time. Since 
they are infrequent and the distribution at small T x — T 
is stationary, one can use the method of characteristics 
to determine the exponentially low probability of forma- 
tion of such a nucleus (instanton for the Fokker-Planck 
equation). An important difference from the theory of 
nucleation in the first-order phase transition is that the 
probability of instanton formation in this case is deter- 
mined by the cooling process. 



III. A TOY MODEL 

To clarify the situation, let us consider the instanton 
solution in the case of one degree of freedom, for which 
the Fokker-Planck equation has the form 



™ = d_ ( D 9W _ vW 
dt dx \ dx 



(13) 



where D is the constant diffusion coefficient and v(x) is 
the macroscopic variation rate of the quantity x with al- 
lowance for its relaxation upon the deviation from equi- 
librium and for an external effect (analogue of phonon 
emission). Setting W — e s , and assuming that the mod- 
uli of S and its first derivative are large, we obtain, to 
leading terms, the equation 



21 -d(—Y- —- — 

dt \ dx J dx dx 



(14) 



This is the Hamilton- Jacobi equation with the Hamilto- 
nian (dS/dx = p) 

The Hamilton equations are the characteristics of this 
equation in partial derivatives, 

dx 



dt 



= -2Dp + v, 



(16) 



dp 
~dt 



dv 
dx 



1>- 



d 2 x 



(17) 



The contribution of the velocity divergence to the Hamil- 
tonian is significant only in the vicinity of the point 
v = 0. We are interested in the special solution that 
passes through the equilibrium point p = 0, v = 0. In 
the ID Fokker-Planck equation, one can eliminate the 
term with a first derivative by substitution; in this case, 
we have an analogy with quantum mechanics and can use 
the well-known results. Nevertheless, we will use direct 
estimates in the vicinity of v = 0. 

In the Hamilton equation, the energy is conserved. In 
view of the smallness of the divergence term, this gives 
H = —Dp 2 + pv = 0, whence p = v/D and 



D 



dx. 



(18) 



We assume that the velocity v(x) is a convex-down func- 
tion with two zeros (stable at zero and unstable at x* 
(x* > 0). Such a shape of the function v(x) is ensured 
by the entire cooling process, including phonon emission. 
For x > x* , the solution tends to larger values of x, while 
the action is gathered from zero to x*, where v < 0. In 
the vicinity of x* , we must take into account the quantity 
dv/dx . For large values of x, p 2 can be ignored, yielding 
p w — (dv/dx) /v, 



S ~ So - ln(v/uo), 



(19) 



where vo is the effective velocity in the region where the 
solutions for x < x* and x > x* match. The solution So 
itself has the form e s sa voe Sa /v, and current j w Voe So . 
One can estimate the value of vo, assuming that all terms 
in the Hamiltonian H are of the same order of magnitude: 



Dp <~ vp 



dv Vr, 



dx x* 



which gives 



(20) 



(21) 



In the many-dimensional case, the situation is the 
same, 




dr l 



dpi 
dt 



dv k d(divv) 
-Pk - 



dx % 



dx 1 



(22) 



(23) 



where p = at the beginning and p — > at the end of the 
trajectory. Consequently, \p\ reaches its maximal value 



4 



somewhere on the trajectory. At this point, the matrix 
dvi/dxk has one zero eigenvalue and p is tangent to the 
corresponding eigenvector; subsequently, the trajectory 
passes to the neighborhood of the point corresponding to 
zero velocity v. This leads to the definition of the critical 
fluctuation (instanton) as a solution passing through the 
point x = p = 0, whereupon p — > for |a;| — > oo as 1/v, 
retaining the probability flux at a constant level. 



IV. OPTIMAL FLUCTUATION 

An analogous procedure can be carried out for the field 
as well. In this case, the Hamiltonian has the form, in 
accordance with Ea. (|lll) . 



H- 



p(r) 



X T o. 



-V 2 p(r) + xV 2 T+U(T c -T) 



T-T r 



Tph 



with the Hamilton equations 
dT 2 X T C 2 



at 



V 2 p{r) + X V 2 T + U{T C -T) 



T-T c 

Tph 



^ = - x V 2 p-^U(T c -T). 
at T ph 



(24) 
(25) 

(26) 



Here, p = 5S/ST(r). Eqs. JjEJ, JU define the critical 
fluctuation and can be reduced to dimensionless variables 
by the substitutions £ = r/^jx T ph, t — t/r p h, 

_ T -T c _ nCpiTpo - T C )U 

rp rp 1 P rp2 ' 

± OO ± C ± OO 

where and n are the new dimensionless fields. In this 
case, the dimensionless equations have the form 



50 

a7 



v 2 e + qu(-b) + 2v 2 n, 



^ = -V 2 n-nu(-e). 

OT 



(27) 



(28) 



The solution should fulfill the conditions II^oo — > 0, 
II r ^-oo — > 0, O r ^_oo — > 1, 0£-»oo — > 1 and pass through 
the neighborhood of $0/<9r « 0, II « at r — > r*. Later 
the fluctuation is developed by cooling, while the random 
fluxes can be neglected and dQ/dr sa V 2 + QU(— 0). 
With exponential precision we can assume that dQ/dr — > 
at r — > +oo, and O tends to the stationary solution st 
of the thermal diffusion equation, Ea. (|2*7|) : 



Q st = sin(0/^, C < 7T , 

6 st = 1 - 7T/e, £ > 7T . 



(29) 



The difficulties with the numerical solution of this 
boundary value problem are due to the instability of 
Eq.l|28p for ascending time whereas Ea. (|27|) is unstable 



for descending time. Therefore, it is impossible to find 
numerically the solution of the Cauchy problem in cither 
direction of time. We briefly describe our adopted pro- 
cedure. At early stages of the evolution, when > 
everywhere in space, it is easy to check the validity of 
the relation 



© = i - n, 



(30) 



which coincides with the thermodynamical theory of tem- 
perature fluctuations. The function II grows with time 
according to Eq. (|28|l (thermal diffusion equation with 
negative time derivative). We can assume that at t = 
the maximum of II will reach 1. After this Eq. (|30() is no 
longer valid. We can consider Eq. (|28|l as a Schroedinger 
equation with imaginary time and II — > at infinite time. 
At large r, the function will be close to the station- 
ary solution, Ea. l|29|l . which becomes zero at £ = ir. This 
means that II, at large r, has an asymptotic proportional 
to 



n m/ = exp(-|A|r)*A, 



(31) 



where ty\ corresponds to the eigenfunction with the 
negative eigenvalue, A = —0.4576, according to the 
Schroedinger equation with the potential —U(tt — £) (all 
other states will grow with r) . 

Let us denote by r* (r) the space-point for which = 
at time r. The curve r* (r) is a function of t starting at 
small t as square root of r (because has a minimum at 
£ = 0) and tending exponentially to ir at large r (accord- 
ing to the Schroedinger-equation analogy, Eq.ffiip). We 
don't know the exact form of r*(r) but we can choose 
some probe function with the same asymptotic behavior. 
Having such a probe function, we can numerically solve 
Eq.(J2BJl integrating it backward in time (it is unstable 
while integrating it forward in time) with the condition 
Hm/ = (of form Eq. i|3*T)l) (the prefactor a should 
be chosen to satisfy Eq. I|30|l at r = 0). 

Using this function n, we can numerically solve Eq. 
(|27|l in two space-time regions independently. The first 
one is the internal region £ < r*(r) the other one is 
£ > r*(r) (the external region). — * 1 for £ — * oo 
(for the external region) and 0(r*) = (for the both 
regions). For the exact function r* the space derivative 
= <90/9£ will be continuous. For our probe function 
r* there will be some jump of the space derivative at 
r*(r). Thus, after performing the calculations, we will 
have in general a nonzero jump-function 

A0' =e'(r"-0)-e'(r*+0) 

depending on the choice of the curve r* . Afterwords we 
should search for a more exact r* in order to minimize 
maa;|A0 |. Some steps within the framework of this pro- 
cedure have been performed. 

We use a space-time grid 2000x2000. The grid spacing 
was taken as <5£ = 0.015 for the space coordinate and 
5t = 0.01 for the time coordinate. First, we perform 
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FIG. 1: Left and right derivatives at r* for two free parameters 



calculations backward in time for II using a standard im- 
plicit numerical scheme where the space derivatives are 
calculated for the final time of each time step. There 
are some modifications of the space grid in the vicinity 
of r = r*(t), for a better finite difference representation 
of the Laplacian. Then we use the analogous scheme for 
Ea. (|27(l and find O in the internal and external regions 
integrating forward in time and obtain AO . We have 
defined our probe function r* (t) by a number of param- 
eters. To obtain an exact solution we need, of course, 
an infinite number of parameters. In practice, for a rea- 
sonable precision we only need a few. The simplest form 
which obeys the asymptotic behavior is 

r* = ay/r + [3t, t < r , 

r* — it — 5exp (— |A|r), r > r . 



FIG. 2: Left and right derivatives at r* for three free param- 
eters 
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FIG. 3: Jump-function AO (r*) for two and three free pa- 
rameters 



The parameters should be chosen such that r* is contin- 
uous and smooth at r = tq. This means that we have 
two free parameters in this case (e.g. a and To). After 
minimization of max\AQ | with respect to these two pa- 
rameters we find right and left derivatives © (r* ± 0) (see 
Fiq.l). 

We can improve our results adding a new term 7T 2 : 



r* = u\fr + [3t + 7T 2 , t < r , 
r* = 7T — S cxp (— |A|t), r > tq. 



(32) 



We have found that for a = 2.28, 7 = 0.016, t = 2.55 
there is an acceptable minimum of | AO | (see Fig. 2). The 
corresponding jump- functions AO for the two and three 
parameter cases (for comparison) are plotted in Fig. 3. 
One can see that the inclusion of this additional term 
decreases the maximum deviation, maa:[|AO |], by a fac- 
tor of 3. 



V. THE OPTIMAL FLUCTUATION 
PROBABILITY 

The solution of Eqs. (|27H . (|28|) allows us to calculate 
the action 

S = / P ^t d3rdt ~ J Hdt 

riCpiT^ — T c ) 3 

= so orp2 (VrvO ■ (33) 

The negative constant sq is the dimensionless action 
00 00 
s = J dr J a r ©nd 3 C - J Htdr, (34) 

—00 —00 

where the space integral is taken over the whole space. 
Here 

Hx= J n(r) [V 2 n(r) + V 2 + J7(-O)0] d 3 £ (35) 
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is the dimensionless Hamiltonian, Eg.l|24|l. The nega- 
tive constant Sq is a universal number corresponding to 
the largest action So and is independent of the values of 
physical constants and the difference — T c . 

From Ea. (|27|l and from the expression for the Hamil- 
tonian, Eq. (J2SJ, it is easy to see that 



dr / nV 2 nd 3 r 



so 



for a smooth solution of Eas. (|27(l . I|28(l . In our case, the 
jump of derivatives should be taken into account which 
modifies the action : 



.s = / dr I nV 2 nd 3 £- I dr I I AG UdS 



can use the relation v = n/T^ . Thus, the number of 
critical fluctuations per unit time in a unit volume is 



3 1/2 

dN n{n(^/xr^K) ) 

dt TpfrCp 



exp 



so- 



nCp (Tf-x 



2T 2 r 



iVXTph) 



The nuclei of a new phase are intensively formed as 
approaches T c and then grow rapidly. We have consid- 
ered the initial phase of critical-fluctuation growth and 
restricted our analysis to the heat transfer via heat con- 
duction, disregarding superfluidity effects at this stage. 
This approximation can be justified by the fact that the 
largest contribution to the action comes from the re- 
gion lying far from the region of low velocities v , where 
p(r) = 5S/5T(r) becomes small and the fluctuation con- 
tribution can be neglected because W ~ v~ 1 e Sa . 



The surface integral should be taken over a sphere of 
radius r* . Using Eq. , 

v 2 n = -<9 T n - u(r* - on, 

we finally find for the action 



s = - / dr I U 2 d 3 £ ~ dt J A0 IldS r « (36) 

— oo £<r* — oo 

Substituting the results of the numerical calculations we 
get 

s = -100.73. 

Test results for a smaller grid spacing, 5£ = 0.003, give 
s = -100.23 

with a negligible jump- function, A0 . These results are 
in agreement with the naive estimate from Eq. l|36|l : The 
characteristic scale of II 2 is 1 and the action is propor- 
tional to the volume of a sphere of radius 7r, which is 
about 100. 

To estimate the temperature-variation rate, one can 
take 



\Vmax \ = — — {WXTphfn. 

Tph 



(37) 



In this case, in accordance with Eq. I)21|l . one can write 
for the probability flux in the transition region 



3,2 



3 



To {n(^/W P ~h) ) 



-e s v. 



(38) 



The constant v cannot be estimated from the theory of 
hydrodynamic fluctuations |l5|. This quantity gives the 
number of small equilibrium fluctuations with ST <C T 
in a unit volume on the atomic scale. As an estimate, we 



VI. BOSE CONDENSATION IN THE 
TWO-DIMENSIONAL CASE 

Most experiments on Bose-condensation in exciton sys- 
tems were done in the two-dimensional (2D) case because 
these excitons are more stable 0,0. Condensation into 
the superfluid state of interwell exctions in AlAs/GaAs 
structures was supposedly observed in In these ex- 
periments excitons were obtained in a 2D quantum well 
and were cooled by phonon emission into the 3D volume 
of the surrounding semiconductor |l8j |. We can expect 
our theory to be suitable to explain the formation of a 
condensate in this case. In such systems the Bose-gas 
of excitons is dilute (nao 2 <ti 1) and has a long lifetime. 
There is no a true condensate at any nonzero tempera- 
ture but it was shown in [nil ] that there is a superfluid 
transition at 



2imh 2 



77ilnln(l/nao 2 ) 



According to [ijj the pressure has the form: 
27m 2 mT 2 C(2) 



m ln(l/nao 2 



2nh 2 



(40) 



(41) 



As it was pointed out in 20] these results hold only under 
the condition, 



lnln(l/na 2 ) > 1. 



(42) 



Using the analogy to the 3D case, we consider tempera- 
ture fluctuations at constant pressure and generalize the 
results of previous sections to the 2D case. In these fluc- 
tuations below T c there is also an increase of the recip- 
rocal phonon time, and in the 2D case Eq.© should be 
replaced by 



Sn _ _ST \n(l/(na 2 )) 
~ ~ ~ ~T {lnln{l/{nao 2 ))) 2 



(43) 
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The large parameter (ln(l/(na 2 ))/(lnln(l/(nao 2 ))) 
plays the role of I/77 3 in Eq. JBJ. Eventually, we ar- 
rive to the same set of equations, Eqs. H27I28J) . as in the 
3D case. Thus, the number of critical fluctuations per 
unit time in a unit volume in the 2D case is 



dN 
~dt 



n(nxT p h) X/2 



exp 



so - 



2T 2 r 



(XTph) 



(44) 

Numerical calculations in the 2D case are somewhat 
more complicated in comparison to the 3D case. There 
is no a stationary solution like Ea. (|29|l . However, if we 
consider the system of a finite size R, there will be an 
analog to Eci. (|2*§)l : 



MO 



Ai ^i(A 1 o)ln(^) 



,. N 7 £ < A 0, 



e. 



(45) 



Here J n is a Bessel function and A x o is the first zero of 
Jo (a;). We see that the negative values of 9 are logarith- 
mically suppressed. This is not very important because 
the system after passing in the vicinity of the stationary 
point evolves further due to the divergence term (which 
was neg lected in Eo.[|2"I)l and Eqs.J23,lEHlO- This term 
in Eq. l(2~4l has the form 



AH 



U(Tc - T) 2 
— a r. 

Tph 



The Hamilton equations, Eas. (|27ll . (|28|l . should also be 
modified: 
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V 2 9 + QU{-Q) + 2V 2 n, 



(46) 



^ = - v 2 n-nc/(-9) 

OT 



s(-e) 



(Too - T c ) nCpXTph 



(47) 



Here 5(— 9) is a delta-function. As in the toy-model, this 
term starts to play an important role at large r while at 
small and intermediate r it is unimportant. There is no 
need to take this term into account in our calculations 
of the action in the 3D case because it is a higher-order 
quasiclassical correction for the 3D instanton. But in 
the 2D case, we can consider the finite-time evolution 
due to this term. The large-time cutoff can be estimated 
by assuming that at this moment the divergence term 
becomes of the same order of magnitude as the main 
terms in Eq. l|47|l . We know the large-time asymptotics 
of n in a finite-size system of radius R, 

a* A (r)exp(-|A|r), 



where ^\(r) is the eigenfunction with the negative eigen- 
value A of the 2D Schrocdinger equation with the poten- 
tial —U(X 1 o — 0- Using this asym 

ptotics and comparing terms in the equation for the 
norm, ||n|| = J n 2 d 2 r, 



Or 
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we can estimate the large-time cutoff as: 
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A 
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The space position of the thermal front, that corresponds 
to this time, is 

Ar = VAt. 

If we choose the point R, where 9(i?) = 1, at the distance 
Ar from the point r* (where 9(r*) = at the largest 
time) , we can assume that the value of the total action 
will be almost independent of the exact position of R. We 
have performed a number of runs with different values of 
R and have found no essential difference in the action. 
The resulting constant is 

.so = —13.6 

The same naive estimate as in 3d can be performed. The 
action now is proportional 7r(A 1 o) which is about 13. 



VII. THE LATER STAGE OF INSTANTON 
GROWTH 

Analysis of the subsequent growth of the instanton re- 
quires the solution of hydrodynamic equations for a su- 
perfluid liquid, because a superfluid core appears in the 
developing fluctuation. We will qualitatively consider the 
phenomena that arise in this case. Proceeding from the 
assumption that the value of T p h is high, we assume that 
the motion in this region is quasi-stationary and tuned 
due to the slow cooling by phonons. We will use the 
hydrodynamic equations for a superfluid liquid in the 
vicinity of the transition point in the form proposed by 
Khalatnikov H3: 



dv s 
dt 



V I y + /J + H S 



Op 
di 



+ div(p s v s -(- p„v„) = 0, 



d d 

g- (PsVl + PnVn) + {pnV l n V k n + p s v\v k s + PS lk ) 



1 — — h 1 div(nerv„) 



= 0, 



dt 
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dt 
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(v„ - v s 
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(52) 



Pa- 



Here, a is the entropy per particle, n is the number of 
particles per unit volume, and p is the density. The sub- 
scripts n, s correspond to the normal and superfluid com- 
ponents, respectively; the constant A is the relaxation 
parameter; and we introduced the term that accounts 
for the phonon-induced energy removal in the equation 
for entropy. Here, the specific chemical potential p s 
for the superfluid density should ensure the condensate 
equilibrium density that is obtained by equating to zero 
the relaxation right-hand side of the equation for p s . In 
our model of a weakly nonideal Bose gas, we can define 
phcnomenologically 

/«. = -^[(»-»(T))] + (48) 

so that p s — m(n — n(T)) — m5n in the equilibrium. 
Here, n(T) is the number of particles outside the con- 
densate. We assume that the quantity Aro/fi is large 
and T c — T is large enough for the approximate equality 
p s + v 2 /2 sa to be satisfied (we disregard quantity v n , 
which is small compared to w s ); this gives 



By introducing the dimensional distance £ — c p r / aur p h 
and v — v s /u we arrive at the equation 



1 -3v 2 



(53) 



The singular points of this differential equation are 

2 1 

£ = 0, v = and £ = -,= , v = -y=, 

the latter point being a focus with the eigenvalues A = 
Since the velocity v must vanish at £ = 0, v « 
£/3 for small £ and increases faster than by the linear 
law, with the derivative || vanishing at v = ^= and 
at a certain £ = £ c , whereupon the derivatives assume 
negative values upon the further increase in v. Thus, 
a regular superfluid flow cannot be continued after the 
point £ c (the constant is on the order of unity and can 
be determined numerically). The critical radius 
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2h 2 a Sn ) ' can ^ e smaner than ^X T ph j it should also be noted 



In this case, it follows from the hydrodynamic equations 
that p ~ p(P, T) = const, 



rp (dna) CpTps 
1 — — h i aiv(nc7v„) = 



(50) 



dt mTph 
and the momentum conservation law gives 
d d 

-g t (psvl + p n vl) + {p n v l n v k n + p s v\v k : + PS lk ) = 

In view of the smallness of v s compared to the sound 
velocity and the smallness of p s , we will neglect these 
corrections to pressure P ~ Pq- In this case, only the 
equation for entropy is significant. Assuming that the 
derivative (dna)/ 'dt is small, according to the assumption 
that the process is quasi-stationary (low temperature- 
variation rate) , we find that the stationary regime 
— div(crnv n ) = —c p ps/T p h should approximately take 
place and that the mass flux should be zero (p s v s + 
Pn.v n = 0). Considering that p n sa p we obtain the equa- 
tion 

CpPs 



-adiv(p s Vs 



Tph 



(51) 



where a is the entropy per particle. This equation de- 
termines the heat transfer in the fluctuation superfluid 
core. Using Ea. 1)49(1 we obtain 



a--r 2 (l 

r 2 dr 



v 2 % c K 



rK -(!--!) 



W Tp h 



0, 



that 1 — v 2 > 0(i.e., a singularity appears in the su- 
perfluid core). This singularity indicates that the quasi- 
stationarity conditions are violated at £ > £ c , and a com- 
plex nonstationary superfluid flow with the intense vortex 
formation in an instanton should appear upon the tran- 
sition to the normal liquid at T > T c . Similar effects are 
observed in a superfluid liquid in the gravitational field, 
where T c is a function of one (vertical) coordinate and a 
fixed heat flux from the superfluid to the normal liquid 
takes place 22]. We are dealing with a similar situation 
arising due to the nonuniform cooling as the critical tem- 
perature in the superfluid nucleus is approached. The re- 
sults of numerical calculation [23j and experimental data 
[24j], [25| indicate the formation of a "vortex" superfluid 
phase with a higher but finite thermal conductivity with- 
out a superfluid transport. The mechanism of vortex 
formation and the vortex phase of this kind have been 
poorly studied both theoretically and experimentally. 



VIII. CONCLUSION. 

Thus, we have shown that, in contrast to 0, a tran- 
sition to the superfluid phase can occur through an in- 
dependent growth of critical fluctuations (instantons) at 
temperatures above the critical point (T > T c ) imme- 
diately in the course of external cooling. These fluctu- 
ations subsequently transform into macroscopic forma- 
tions. The growth of the nucleus of the superfluid state 



9 



is accompanied by vortex generation in its external part. 
Consequently, vortex defects appear both due to the in- 
dependent nucleation with an arbitrary phase upon cool- 
ing (the Zeldovich-Kibble hypothesis) and directly dur- 
ing the growth of each superfluid nucleus. This vortex- 
generation mechanism during the growth of an instanton 
significantly differs from the mechanism determined in 
[26j , where the existence of a superfluid flow interac ting 
with the heated normal regions was presumed. In [26j, 
an attempt was made to explain the results of experi- 
ments 27], in which 3 He was irradiated by neutrons. As 
a result, some regions heated to temperatures above T c 
appeared. These regions were cooled by the surround- 
ing superfluid 3 He, and the formation of vortices was de- 
tected. Thus, nonuniform cooling took place that differs 
considerably from the model used in our study. In the 
critical fluctuation considered here, heating takes place 
due to its nonsuperfluid surroundings. Consequently, it is 
advantageous for the fluctuation to preserve its spherical 
symmetry to reduce this heating. In the case of cooling 
of a heated region with superfluid surroundings |27| , the 
interface must obviously be unstable against its shape 



distortions, because this leads to a faster cooling. How- 
ever, the stability, as well as the phase-transition mecha- 
nism itself, under such conditions (which, in contrast to 
|26j , are not associated with the existence of an external 
superfluid flow) calls for detailed investigations. 

Probably, analogous schemes can be developed for the 
kinetics of various other phase transitions in the presence 
of the external cooling with the scenario essentially given 
by the toy model described in the text. 

IX. ACKNOWLEDGMENTS 

We are grateful to V.V. Lebedev, I.V. Kolokolov and 
H. Miiller-Krumbhaar for discussions. This study was 
supported by the President of the Russian Federation 
(grant no. NSh-1715.2003.3) in Support of Young Rus- 
sian Scientists and Leading Scientific Schools, the pro- 
gram Quantum Macrophysics of the Presidium of the 
Russian Academy of Sciences, and RFFI grants no. 03- 
02-16012, 05-02-16553. 



[1] R.Becker, W.Doering. Annalen der Physik 24 , 719 
(1935) 

[2] Ya. B. Zel'dovich, Zh. Exp. Teor. Fiz. 112, 525 (1942) 

[3] J.S. Langer. Ann. of Physics 54 258 (1962) 

[4] I. M. Lifshits, Zh. Exp. Teor. Fiz 42, 1354 (1962) [Sov. 

Phys. JETP 15, 939 (1962)] 
[5] Ya. B. Zel'dovich, I. Yu. Kobzarev, and L. B. Okun', Zh. 

Exp. Teor. Fiz. 67 , 3, (1974) [Sov. Phys. JETP 40, 1 

(1975)] 

[6] T.W Kibble, J.Phys. A9, 1387 (1976) 

[7] W.H. Zurek, " Cosmological experiments in condensed 
matter system", Physics Reports 177-221 (1996) 

[8] E.A. Brener, S.V. Iordanskiy, R.B. Saptsov, JETP Let- 
ters 79, 410 (2004) 

[9] H. Henry and Herbert Levine, PRE 68, 031914 (2003) 
[10] M. Marder, PRE 54, 3442 (1996) 

[11] A.I. Chernykh, M.G. Stepanov, PRE 64, 026306 (2001) 

[12] L. D. Landau and E. M. Lifshitz, Course of Theoretical 
Physics, Vol.5 "Statistical Physics", Fizmatlit, Moscow, 
2002; Pergamon Press, Oxford, 1980. 

[13] E. M. Lifshitz and L. P. Pitaevskii, Course of Theoreti- 
cal Physics, Vol.9 Statistical Physics p. 2. h. 9, Fizmatlit, 
Moscow(2002); Pergamon, New York, 1980 

[14] V. I. Klyatskin, Stochastic Equations by the Eyes of a 
Physicist, Fizmatlit, Moscow (2001) 



[15] E. M. Lifshitz and L. P. Pitaevskii, Physical Kinetics 
(Nauka, Moscow, 1979; Pergamon Press, Oxford, 1981). 
[16] L.V. Butov, A.I. Filin. PRB 58, 1980 (1998) 
[17] A.V. Larionov et. al. JETP Letters 75, 570 (2002) 
[18] S.V. Iordanskiy, A.B. Kashuba. JETP Letters 73, 542- 
545 (2001) 

[19] V.N. Popov , Functional Integrals in Quantum Field The- 
ory and Statistical Physics, Atomizdat, Moscow (1976). 
(Reidel,Dordrecht,1983). 
[20] D.S. Fisher, PC. Hohenberg. PRB 37, 4936 (1988) 
[21] I. M. Khalatnikov, The Theory of Superfluidity (Nauka, 

Moscow, 1971), Chap. 9 
[22] Akira Onuki. J. Low Temp. Physics 50, 5/6, 433 (1982) 
[23] P.B. Weinman, J. Miller. J. Low Temp. Physics 119, 155 
(2000) 

[24] Feng Chuan Lui, Guenter Ahlers. PRL, 76, 8, 1300 
(1996) 

[25] H. Baddar, G. Ahlers, Kuehn, H. Fu. J. Low Temp. 

Physics 119, 1/2, 1 (2000) 
[26] I.S. Aranson, N.B. Kopnin, V.M. Vinokur. PRL 83, 2000 

(1999) 

[27] V.H.M Ruutu, V.B Eltsov, A.J Gill et al., Nature 382, 
334 (1996) 



